Langevin equation for the extended Rayleigh model with an asymmetric bath 



Alexander V. Plyukhin and Jeremy Schofield 
Chemical Physics Theory Group, Department of Chemistry, 
University of Toronto, Toronto, Ontario, Canada M5S 3H6 
(February 2, 2008) 

In this paper a one- dimensional model of two infinite gases separated by a movable heavy pis- 
(N . ton is considered. The non-linear Langevin equation for the motion of the piston is derived from 

' first principles for the case when the thermodynamic parameters and/or the molecular masses of 

gas particles on left and right sides of the piston are different. Microscopic expressions involving 
time correlation functions of the force between bath particles and the piston are obtained for all 
parameters appearing in the non-linear Langevin equation. It is demonstrated that the equation has 
stationary solutions corresponding to directional fluctuation-induced drift in the absence of system- 
atic forces. In the case of ideal gases interacting with the piston via a quadratic repulsive potential, 
the model is exactly solvable and explicit expressions for the kinetic coefficients in the non-linear 
• Langevin equation are derived. The transient solution of the non- linear Langevin equation is ana- 

lyzed perturbatively and it is demonstrated that previously obtained results for systems with the 
hard-wall interaction are recovered. 
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& . I. INTRODUCTION 

The Brownian motion of a massive piston in a cylinder filled with an ideal gas is one of the oldest models of 
non-equilibrium statistical physics. In its simplest version, both the piston (of mass M) and the gas molecules (of 
mass m) are confined to move in one dimension along the symmetry axis of the cylinder. The piston is assumed to 
| be adiabatic in a sense that its heat conductivity is negligible so that there is no transfer of heat between the two 
<—j . compartments of the cylinder as long as the piston is held fixed. This simple model has served as a useful test of new 
schemes and approximations for several generations of physicists (see Ref. [1] for references). In most treatments it is 
assumed that the bath particles interact with the piston via the hard-wall potential, i.e. piston-molecule collisions are 
considered to be instantaneous. In this model the piston interacts with bath particles through a sequence of binary 
7-H \ collisions and the possibility of simultaneous interaction of the piston with more than one molecule is neglected. 
■ This approximation is valid as long as one is interested in the asymptotic long time behavior of the system which is 
apparently insensitive to details of the interactions between the bath particles and the piston. On the other hand, 
the motion of the system on short to intermediate time scales is clearly influenced by the nature of bath particle- 
piston interactions, and one has to take into account the finite range of the potential and effects of multi-particle 
collisions. It is known that many-particle collisions are an important factor in liquid-like systems, essentially affecting 
^•f-) ' the shape of the velocity correlation function, especially when a repulsive part of potential is relatively "soft" [2] . In 
such systems, a given particle interacts simultaneously with a relatively heavy swarm of other particles in its vicinity, 
leading to an enhancement of caging effects in comparison to systems with a short-ranged potential. As a result, 
the velocity autocorrelation function has a more pronounced negative part corresponding to an anti-correlation in 
the velocity induced by a particle rebounding of its neighbors. Multi-particle collisions may also be important for 
tagged particle motion in a dilute gas-like medium provided the tagged particle (or "piston" in the present case) is 
large enough that the average number N of gas molecules in the interaction shell around the particle is larger than 
one. Although in this case the qualitative form of correlation functions does not depend on TV, the softness of the 
potential essentially influences the short-time behavior of the system leading to the non-exponential initial decay of 
time correlation functions [3]. 

. ^ \ Another reason why one may wish to go beyond the approximation of hard- wall interaction is methodological. 
It is often problematic to directly apply general analytical techniques involving differential operators for systems 
with a singular potential. Even if one expects multiple collisions to be unimportant, it is often convenient to model 
interactions via a short-ranged potential of range U and subsequently analyze results in the hard- wall limit — > 0. 

In a earlier paper [1], a microscopic derivation of the non- linear Langevin equation was presented for a system 
consisting of a massive piston interacting with a bath particles via an arbitrary repulsive potential. We use the 
term "non-linear" to refer the Langevin equation with not only the linear dissipative term (Stokes damping) but also 
dissipative terms of higher orders in particle's momentum. In Ref. [1] all generalized kinetic coefficients appearing 
the non-linear Langevin equation were expressed in terms of time correlation functions involving the interaction force 
between the piston and the bath particles and its derivatives. For some systems, such as an ideal gas interacting 
with the piston through a repulsive quadratic potential, the microscopic expressions for all the kinetic coefficients can 
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be calculated analytically. For this reason, simple models, such as the ideal gas-piston system, provide a convenient 
means of studying many subtle points of Brownian motion theory, such as the role and form of non-linear damping, 
relative importance of non-Markovian effects, and convergence properties of small parameter expansions. 

In Ref. [1] the dynamics of the piston in a homogeneous bath was examined. The purpose of this paper is to extend 
the analysis in Ref. [1] to the case of an asymmetric bath in which the parameters characterizing the bath to the left 
and to the right of the piston are different. This model has received a recent renewal of interest after it was discovered 
that the system may exhibit non-trivial transient and stationary behavior when there is an initial asymmetry in the 
thermodynamics properties of the gases to the left and to the right of the piston [4]. Among other interesting points, 
it was found that the piston undergoes a noise-induced directional movement in the absence of macroscopic forces, 
a characteristic of motion in molecular motors and stochastic ratchets. In the limit of an infinitely long cylinder, 
Gruber and Piasecki found [5] a stationary solution of the equations of motion of the piston corresponding to the 
drift of the piston in the direction of the compartment with higher temperature even when the pressure on the left 
and the right of the piston are the same. This result has been later confirmed by numerical simulations [6]. This 
fluctuation-induced motion is an effect of the first order in a small mass-ratio parameter A, defined as A = y/m/M. 
It does not follow from the conventional Langevin equation with linear dissipation, and appears only when non-linear 
damping terms, which are of higher orders in A, are taken into account. 

In this paper the projection-operator method applied by Mazur and Oppcnhcim [7] to the theory of Brownian is 
adapted to the case of the asymmetric bath. We shall consider a slightly more general model of the asymmetric bath 
than that analyzed by Gruber and Piasecki in that not only temperatures but also the masses of the bath particles in 
the left and right compartment may differ. For such a system, the Langevin equation, including non-linear dissipative 
terms to third order in A, is derived from first principles. The resulting equation is general, holds for arbitrary 
interactions between the piston and the bath particles, and is not restricted to the ideal gas bath. The kinetic 
coefficients appearing in the non-linear Langevin equation are expressed as integrals of time correlation functions of 
the force between the piston and bath particles. For a bath of ideal gas particles interacting with the piston via the 
parabolic repulsive potential, the correlation functions can be computed in closed analytical form. For this model, an 
explicit expression describing relaxation of the momentum of the massive piston is obtained. The analysis based on 
the Langevin equation is much more simple than that involving the language of distribution functions adopted in other 
papers on the subject. The perturbative and stationary solutions are analyzed and demonstrated to be consistent 
with the results of Gruber and Piasecki obtained for a model of instantaneous binary collisions. 



II. AN EXACT EQUATION OF MOTION FOR THE PISTON 

The system consists of a piston of mass M and cross-sectional area S confined to move in one dimension (chosen 
to be along the x axis) in a cylinder of total length L. It is assumed that the piston is initially at position X near 
the center of the cylinder, taken to be the origin, with the leftmost and rightmost-ends of the cylinder at positions 
—L/2 and L/2, respectively. The left and right compartments of the cylinder arc filled with bath molecules of density 
ni and n r , temperature Ti and T r , and mass mi and m r , and move in one dimension along the x axis. We assume 
the bath molecules are confined to their respective compartments by hard-wall interactions with the immobile ends 
of the cylinder, and interact with the piston through a short-ranged, repulsive potential. The Hamiltonian of the 
piston-bath system can be written in the form 

H=^ + H (X), (1) 

where X and P are the coordinate and the moment of the piston, and Ho(X) is the Hamiltonian of the bath 
in the presence of the piston constrained at position X. H (X) can be decomposed into the sum of two parts, 
H a (X) = Hq(X) + Hq(X), corresponding to the Hamiltonian for the bath molecules on the left and right of the 
piston, 



1=1 K J 



(2) 



For simplicity of notation, above and throughout this paper the super- or subscript index a. = {I, r} is used to label 
dynamical variables in the left and right compartments of the cylinder. In the above equation N a and m a are the 
number and the mass of particle in a respective compartment. 

In the following, we restrict our analysis to the case when differences in the thermodynamic parameters of the bath 
in the left and right compartments of the cylinder are small. Furthermore, the molecular masses mi and m r are 
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assumed to be of the same order of magnitude and much less than the mass of the piston M. Under these conditions, 
one might anticipate that the directional contribution to the momentum of the piston is small, and on average is of 
order P ~ y/MksT p , where the effective temperature T p of the piston is of the same order of magnitude as the left 
and right temperatures 7} and T r . In the subsequent analysis, we show that this intuition is correct, and an explicit 
expression for T p will be presented. 

It is convenient to express the equa tions o f motion of the piston-bath system in scaled coordinates. To this end, we 
introduce the small parameter A = y/m/M, where m is an arbitrary mass of the same order of magnitude as m; and 
m r . One may reasonably expect that the scaled momentum of the piston P* = XP will typically be of the same order 
of magnitude as the momentum of a bath molecule. The parameter A therefore serves as a quantitative measure of 
the time scale separation between the slow evolution of the massive piston and the fast evolution of the light bath 
molecules. As in the case of Brownian motion [7], it will be useful as an expansion parameter to simplify the physics 
of system [8]. 

In terms of the scaled piston's momentum P* = XP, the Hamiltonian reads H = P^ /2m + Ho. The corresponding 
Liouville operator C, which governs the evolution of an arbitrary dynamical variable via the equation A{t) = e ct A 
(A = A(0) throughout the text), can be written in the form 



Here the Liouville operators Cq 



C — Cq + XL\, £q — £q + £q. (3) 



describe dynamics of the bath molecules in the left and right compartments of the cylinder in the presence of the 
piston held fixed at position X, where the force Fi acting on a molecule i is F, = —dU(X — Xi)/dxi. The Liouville 
operator C\ is defined by 

where F = — ^\ dU(X — Xi)/dX is the total force exerted on the piston by the bath molecules. 

Our task is to express the equation of motion of the scaled momentum of the piston, dP*(t)/dt = XF(t), in the 
form of a Langevin equation in which the instantaneous force F(t) acting on the piston at time t is decomposed into 
"dissipativc" and "random" parts. To accomplish this goal, we define the projection operators Pi, V r and V, which 
act on an arbitrary dynamical variable A according to the rules 

P a A= (A) a = J dQ aPa A, (6) 

PA = P{P r A = (A) = J dQt pi j dfl r Pr A, (7) 

where Q = {xf,pf} are the phase points for the bath molecules in the left and right compartments of the cylinder, 
pa = Z~ 1 exp{—f3 a Ho(X)} are the corresponding canonical distributions for the bath molecules in the presence of the 
piston fixed at coordinate X, j3 a = (fc^Ta) -1 an d fcs is the Boltzmann's constant. Note that the projection operator 
P effectively averages over initial conditions of the bath at a fixed position of the piston. Using the operator identity 

e (A+B)t = £ At + r* dre A(t-T) Be (A+B)T^ (g) 
JO 

with A = C and B = —PC, one obtains the following decomposition for the force on the piston F(t) — e ct F, 

F(t) = Ft(t) + / dr e^-^PCF^r), (9) 
Jo 

where F*(t) = e QCt F and Q = 1 - P. 

The structure PCF^(t) in Eq.(9) can be further simplified to XP C\F\t) due to the orthogonality property, 

PtC l - P r Cl = PCo = 0, (10) 
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which follows from the relation C%p a = 0. Using this relation and the definition of the projection operators, Eq. (9) 
takes the form 

F( t ) . F>(t) + A [ * + ^(FF'M) } . (11) 

In this expression the factor 

(Jx F ^ T )) = J dfl lPl J dn rPr ^FHr), (12) 
can be simplified by pulling the differential operator out of the integral, 

(^ t(T) ) = M^^'j dQl J d ^F\r)^(p lPr ). (13) 

Using properties of the projection operator, the first term on the right-hand side of Eq. (13) can be simplified using 
the fact that 

(F\t)) = (e tQC F) = Ve tQC F = VT = (T) . (14) 

For a homogeneous system, one expects (F) = and hence (F^(t)) = 0. For an asymmetric system, on the other 
hand, (F) ^ and the first term on the right-hand side of Eq. (13) generally does not vanish. 

The second term in the right-hand side of Eq. (13), involving -J^pip r , can be worked out taking into account 
that in the distributions p a = Z~ 1 e~ l3aH « not only the Hamiltonians but also the partition functions Z a depend 
parametrically on X, 

d ( 1 QZ 1 dZ \ 

gx(piPr) = PWr [fkFi + (3 r F r - (15) 



It is straightforward to show that 

1 dZ a 



Z a dX 



= P a (F a ), (16) 



where Fi and F r are the forces on the left and right surfaces of the fixed piston such that Fi + F r = F. Therefore, 
Eq. (13) takes the form 

(^ FHt) ) = ^-MFiFHm-MFrFHt))), (17) 

where the double brackets denote the cumulants, ((AB)) — (AB) — (A)(B). Substituting this result into Eq. (11), we 
get the following exact equation of motion of the piston 

^ = XFHt) + A 2 f dr e^) l^FF^r)) + ^ 
dt w J \dP* V " m dX 

~m^{r))) - Ml((F r Ft (r )))J . (18) 

For some systems, such as the extended Rayleigh model discussed in Section 4, the parametric dependence of (F) 
on the position X of the piston is weak, and the derivative of (F) scales as 1/L in the limit of large L. In fact 
for the Rayleigh model, the bath is comprised of ideal gas molecules, and the parametric dependence of (F) on X 
arises through the dependence of (F) on the concentrations of bath molecules in the respective compartments of the 
cylinder, n a . Since 

Ni N r 

n l = TTTT, TTTFi n r 



(L/2 + X)S r {L/2-X)S' 

where Ni and N r are the total number of bath molecules in the left and right compartments, respectively, dn a /dX ~ 
n a /L for large L, and hence d(F)/dX ~ 1/L. Thus, for large cylinders, this term may be dropped in the equation of 
motion of the piston and we obtain 
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= Aft(t) + A 2 fdr e A*-) J ^ (FF t (r)) _M ((F;F t (r))) 



dt 



which is exact for such systems in the thermodynamic limit. By expanding Eq. (19) in powers of the square root of 
the mass ration A, one can derive the Langevin equation to any order in A. 



ML((F r F\r)))\ , (19) 



III. THE NON-LINEAR LANGEVIN EQUATION 

The force F^(t) = e QCt F in the Eq.(19) can be expanded in powers of A using the fact that V£o = and the 
operator identity (8), to yield 



F f (i) = e^a+XQC^t F 



= F (t) + A f dre^-^QdFoir) + 0(A 2 ), 
Jo 



(20) 



where F a (t) = e Cot F is the force exerted by the bath molecules on the fixed piston. Since CiF (t) = "^^-^oW, the 
above equation takes the form 

A P /"* 

F*(t) = Fo(t) + — dT5G(t,r)+0(X 2 ), (21) 
m Jo 

where 

G(t,T) = e^ t -T)jL Fo (T) (22) 

and SG = G — (G). Using this result, one can extract the leading order approximation in A of the correlation 
functions (FF^(t)} and ((F a F^(t)}} appearing in the equation of motion (19). Since the evolution of the dynamical 
functions F (t) and G(t) is governed by the Liouville operator £ , the left and right components of F (t) and G{t) 
are independent, which implies that the cumulants ((FiFq)}, ((F r Fo)} 7 ((FiG r }} 7 and {{F r Gi}} are all zero. Taking this 
into account, one obtains 

((F a F\t))) = ((F a F a (t))) + —^ j dt'((F a G a (t,t'))), 

and 

(FF^t)) = (FF Q (t)) + ^* f dt'lmG^t 1 ))) + ((F r G r (t,t')))\. 

m J Q I > 

Substitution of these expressions into the equation of motion (19) leads to the non-Markovian (generalized) non-linear 
Langevin equation, 

d ^fi = XFHt) - A 2 [ drP^t - r)Mi(r) - A 3 /* drP?(t - r)M 2 (r) 
dt Jo Jo 

+A 3 f dTM 3 (r) +0(A 4 ), (23) 
Jo 

where the memory functions are given by 

Ml(T) = ^{MFi*1(t)))+MFtF5(t)))}, (24) 
M 2 (r) = ^r / dT'{(3i((FiGi(r,r')))+0 r ((FrG r (T,r')))}, (25) 

M 3 (t) = - f dr'hF^T'))) + «F r G r (r, /)))}. (26) 
m Jo I > 
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Note that for a totally symmetric system in which all parameters characterizing the bath to the left and to the right 
of the piston are the same, ((FiGi(t,t'))) = — ((F r G r (r, r'))). In this case the functions M 2 (i), M 3 (t) vanish, and 
the Langevin equation (23) is linear. Furthermore, for a symmetric bath it can be shown that the first non-linear 
correction term is of order A 4 and proportional to P 3 [1]. 

Assuming that the memory functions Mi (t) decay with a characteristic time r c which is short on the time-scale for 
relaxation of the momentum of the piston, the generalized Langevin equation (23) can be written in form that is local 
in time. In fact, on a time scale r < r c the momentum changes primarily due to the random force while the effect of 
the damping force is of higher order in A, 



P„(i - r) = P.(t) - f dr'A(r') = P* (t) — A f dr 1 F\t') + 0(\ 2 ). 

Jt-T Jt-T 



(27) 



Substitution of this expression and a similar one for P 2 (t — r) into Eq.(23) leads to an expression that is local in time 
and of the form, 

^> = AFt(t) - A 2 G(t)P*(t) - \%(t)P?(t) + \ 3 ( 3 (t) + G(A 4 ), (28) 
where Q (t) = J * dr Mi (r) and (t) is the random force modified with a correction of second order in A 



PUt) = FUt) - — f drMAr) [ dr'F^T 1 ). 
m Jo Jt-T 



(29) 



Since this correction is small and docs not change statistical properties of the force, it will be neglected below. In 
what follows the dynamics of the piston on a time-scale longer than the bath correlation time t c are considered. When 
( > t Ci the kinetic coefficients Q(t) assume their limiting values Q = J Q °° drMj(r). Re-expressing this equation in 
terms of the unsealed momentum of the piston P = P*/A yields 

^ = Ft(t)- 7l P(t)- 72 P 2 (t) +73 , (30) 
where the kinetic coefficients 7i are given by 

71 = A 2 Ci = ^ f°° dt^mF^r))) + p r ((F r FZ(T)))}, (31) 

72 = A 4 C 2 = ^p J dtj dr^dFlGlit, t)» + Pr((F r G r (t, r)))}, (32) 

73 - A 2 C 3 = ± J™ dt J ^{((FtGtfrT))) + ((F r G r (t,r)))}. (33) 

It is important to note that Eq. (30) with the kinetic coefficients given by Eqs. (31) - (33) is a general result that 
is valid for arbitrary interaction potentials for the bath and that is also independent of the specific form of the 
interactions between the piston and the bath molecules. 

Note that the correlations ((FiGi(t, t))) and ((F r G r (t, r))) are of different sign, which suggests that for the asymmetric 
bath either 7 2 or 7 3 can vanish for certain combinations of the parameters of the bath. If 7 2 = 0, the Langevin equation 
is linear and can be easily integrated. We postpone further discussion of these issues until the next section where 
the general results (30)- (33) are applied for the specific model of ideal gas molecules interacting with the piston via 
a truncated parabolic potential. 



IV. THE EXTENDED RAYLEIGH MODEL 



The kinetic constants 7 j in Eqs. (31)-(33) are expressed in terms of integrals of correlation functions of dynamical 
variables whose time evolution is described by the Liouvillian operator Co of the bath in the presence of the fixed 
piston. These correlations can be calculated explicitly [1] in the case when the bath is comprised of ideal gas molecules 
interacting with the piston via a parabolic repulsive potential 

r \k f ( x - xtf , x >xt 

U{x-X)=l \k f {x-X r f , x<X r (34) 
[ , otherwise, 



G 



Here x is the coordinate of a bath molecule, Xi = X — a, X r = X + a define the boundaries of interaction between 
the molecule and the piston, with the parameter a serving as a measure of the range of the potential. In the model 
system, the width of the piston is neglected. In addition, it is assumed that the temperatures of the ideal gas in 
both compartments arc sufficiently low (or the potential strength constant kf is sufficiently large) that the average 
penetration (fcj/3 a ) -1 / 2 of a bath particle into the interaction regions [X a ,X) is much less than a [9]. 

In the limit of large cylinder length L, rccollisions of the piston and gas due to the finite size of the bath can be 
ignored. This assumption corresponds to analyzing the motion of the piston on intermediate time scales t <C tl, 
where tl is the time it takes on average for a bath particle to travel half the length of the cylinder. For t <C tl, the 
force acting on the sides of the piston can be written in the form [1] 

Fi(t) = -k f [ dv f dqN(Xi + q,v;t-n)— sin^, 

J0 J-VT, v 

f° f~ VTr v Ld r q 

F r (t) = —kf / dv dqN(X r + q,v;t — r r ) — sin— — , 

J-oo JO w r v 



where oj a = y/kf/m a is the characteristic frequency of the parabolic system, r a = ir/cu a , and the function N(x,v;t) 
is the microscopic linear density of particles defined by 

N(x, v;t)=J2 S(x - Xi(t))S(v - Vi(t)). 

i 

The average forces are: 

<**> = ¥ <^>=-¥- (35) 

PI Pr 

Note that as mentioned in Section 2, the average force (F) = (Fi) + (F r ) does not vanish in general in the asymmetric 
system though the derivative d(F) / dX ~ 1/L can be neglected for a long cylinder. The equation of motion of the 
piston is given by Eq. (19) in this limit. 
For this model, it was found that [1] 

((F a F a (t))) = Jkf~/i a \(F a )\ (t), (36) 
{{F a G a (hM))) - -k f (F a ) &{hM), (37) 

where the functions £,f(t) are 

£i(t) — —^=lsmu a t + (it — uj a t) cosLd a t\ 9(r a — t), (38) 



2tt 

(*i)*2) = ^ cosw Q i 2 (l + cosw Q ii)0(r a - t x ) 9(T a - t 2 ), (39) 

and 6(t) is the step function. 

The above results and Eqs. (31)-(33) lead to the following expressions for the kinetic coefficients, 

(40) 



72 = -j^{mifl(F|)+mrW)}, (41) 
^{ mi ( Fl )+m r (F r )}. (42) 



73 = M 

Note that these expressions are independent of the force constant kf of the quadratic potential. 

Since (F a ) = ±n a S/f3 a one can see from (41) that 72 = when the mass densities in both compartments are the 
same, mini = m r n r . In this case the Langevin equation (30) becomes linear and has a solution 



P{t)=Pe-^ t + / ^ 1 (*- T ) J F lt (r) + ^ 3 -(l-e-T lt ) 
Jo 7i 



(43) 
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which describes relaxation of the momentum to the stationary value (P) = (P) /71 + 73 /71 . If (P) = it follows from 
(40) and (42) that the stationary momentum equals 



m ^ mi (44) 

which means that the piston moves in the direction of the compartment with the heavier bath particles. Conditions 
of equal mass densities, m ; n; = m r n r , and of equal pressure, n;//3/ = n r /f3 r , are satisfied simultaneously only when 
rriifii = m r (3 r . Therefore, Eq.(44) describes motion in the direction of the compartment of higher temperature. 

In the general case where 72 7^ 0, the non-linear Langcvin equation (30) has a form of the non-linear Riccati equation 
and cannot be explicitly integrated. Using the transformation u = exp(7 2 P) it can be converted into a second-order 
linear equation. In the following, however, this line of reasoning will not be pursued and a simple perturbation analysis 
of the non-linear Langevin equation will be conducted. 




V. STATIONARY AND TRANSIENT SOLUTIONS 

Let us first analyze the stationary solution of Eq. (30). Since (F^(t)) = (F), as demonstrated in Eq. (14), the 
Langevin equation (30) implies that the stationary value of the momentum of the piston is given by 

(P) = 1<F)_ J(P 2 ) + J. (45) 
7i 7i 7i 

To calculate (P) perturbatively to first order in A, one substitutes in Eq. (45) the lowest in A approximations for (P 2 ), 
which can be derived solving the linear Langevin equation 

P{t)=F {t)- ll P(t). (46) 

Note that here the "random" force F generally is not zero-centered. In the long time limit, the auto-correlation 
functions for the force on the piston in the left and right compartments are effectively 

(F a F»(t)) = (P Q > 2 + ((F a F°(t))) -+ (P Q > 2 + 2r a <J(t), (47) 
where T a = J °° dt ((F a F®(t))). Thus, the correlation function of the total force F = Pq + Pq is 

(PPo(t)) = (F) 2 + 2TS(t), (48) 

where r = T ; +T r . 

Using (48), one can calculate (P 2 (t)) from the solution of the linear Langevin equation (46), namely, 

P (t) = Pe- 7lt + / dTe-> l(t - T) F (T). (49) 
Jo 

For i> 1 /71 , the result assumes the form 

CPo 2 > = - + ^> 2 , (50) 

7i 7i 

where T is given by 

poo />oo 

T= dt((F l F?(t)))+ dt{{F r F?(t))). (51) 
Jo Jo 

It has been assumed throughout our analysis that the contribution of the systematic force (P) to the momentum 
of the piston is small. This implies that the second term on the right hand side of Eq. (50) is small compared to the 
first one. Since r ~ A and 71 <~ A 2 this restriction requires that (P) <~ A 1+c , for some e > 0. 

For the truncated quadratic potential, T can be evaluated using Eq. (36) to obtain 
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It then follows that for the case (F) = 0, one obtains 



(P 2 )=M ^^ + ^ -SS. (53) 

For a random force obeying Gaussian statistics, distribution of the momentum of the piston is of Maxwcllian form 
with the effective temperature fc^Tp = (Pq)/M. In particular, if to; = m r , then it is found that T p — \JT{F r from 
Eq. (53), in agreement with the results of Gruber and Piasecki [5]. 

Substitution of expression (50) for (Pq ) into Eq. (45) yields for the stationary momentum, 

(P) = — (F) - ^| (F) 2 + ^ - (54) 
7i 7i 7i 7? 

For (F) = and with explicit expressions (40)-(42) for ji, Eq. (54) takes the form 

When the masses of the bath particles are the same in both compartments, to; = m r = to, this expression reduces to 
the result of Gruber and Piasecki [5] 




(P) = y|V^(v / M : -V / M). (56) 

Note that in this case, 73 vanishes as can be seen from Eq. (42). If m;/3; = m r (3 r (and therefore 72 = 0), Eq. (55) 
coincides with the result (44) obtained in the previous section. 

In addition to the stationary solutions of the non-linear Langevin equation, perturbative solutions of the time 
evolution of the momentum of the piston can be examined by expanding the momentum in the form P(t) = Po(t) + 
Pi(t) + - ■ ■ where Po(t) is the solution (49) of the linear Langevin equation (46) and Pi(t) is the next-order correction 
in A which satisfies the equation 

= Fl (t) - 7iPi(i) - 72 ^o 2 (*) + 73, (57) 

and the condition Pi(0) = 0. In this equation, the force F\(t) = F^(t) — Fo(t) is zero centered, since (F^(t)) = 
(F (t)) = (F(0)), as established in Eq. (14). Solving Eq. (57) and taking the average gives 

(P^i)) = (1 - e -Ti*) _ 2£ F 2 e - 7l t ^ _ e - 7l t) 

-^l(F)Pe-^ t { ll t-(l^e-^ t )\ 
7i L ' 

-^r(^) 2 e~ 7lt (sinh 7l i - 7l t) 



2 72 r 



7? 



-e 



" 71 * (cosh 7 it- 1) , (58) 



where P = P(0). Note that the average momentum (P) = (Po(t)) + (Pi(t)) takes the stationary value (54) in the 
long-time limit t > 7^. It should be emphasized that the transient solution obtained here is valid only on a time scale 
that is much longer than the characteristic time t c for the relaxation of the bath but shorter than tl , the time scale 
for particles to move half the length of the cylinder. The above analysis may be readily extended to the short-time 
domain using as a starting point the Langevin equation (28) with time-dependent damping coefficients. 
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